Detection of long-lived staus and gravitinos at the ILC 
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Abstract. A study is presented illustrating the excellent potential of future International Linear 
Collider (ILC) experiments to detect metastable staus f , measure precisely their mass and lifetime, 
and to determine the mass of the gravitino G from the decay f — > tG, thus providing direct access 
to the gravitational coupling, respectively Planck scale. 



1 Introduction 

Supersymmetry (Susy) provides an attractive scenario 
to account for the amount of dark matter in the uni- 
verse. If i?-parity is conserved, the lightest supersym- 
mctric particle (LSP) is stable and an ideal dark mat- 
ter candidate. A very interesting option is the spin 3/2 
gravitino G. The mass of the gravitino is set by the 
Susy breaking scale F via m^/ 2 = m G 

= F/V3Mp, 

with M P ~ 2.4- 10 18 GcV the reduced Planck scale. In 
general m 3 / 2 1S a f ree parameter and may extend over 
a wide range of 0(eV — TeV) for gauge, gaugino and 
supcrgravity mediated symmetry breaking. 

A gravitino LSP may be produced in decays of 
Susy particles. If the next-to-lightest supersymmet- 
ric particle (NLSP) is the scalar tau f , the dominant 
process is f — > tG. Since the coupling is gravitational, 
the lifetime may be very long, ranging from seconds to 
years. The decay- width If, respectively lifetime r = 
rr\ of the t NLSP 

* r 2 1 4 
r - = : —L. i Q_ i 

f ^ rG 487rMj,m| m? [ ' 

depends only on the masses rrif and m,Q as well as on 
the Planck scale Mp - no further SUSY parameters are 
required. 

The cosmological production of gravitino dark mat- 
ter proceeds essentially via thermal production and/or 
late decays of the NLSP. The big bang nucleosynthesis 
puts constraints on the f lifetime [1], e.g. from the en- 
ergy release in hadronic decays one expects r < 10 7 s 
for itIq ~ 100 GeV. Bound states of Mf~ may alter 
the production of light elements considerably, although 
possible consequences are controversial. To avoid 6 Li 
overproduction the f lifetime should be restricted to 
r < 5 • 10 3 s [2], while the author of [3] argues that the 
prior synthesised elements may be destroyed again at 
longer lifetimes. 

Experiments at the ILC offer a unique possibility 
to detect long-lived staus and to study the proper- 
ties of gravitinos, which otherwise cannot be observed 



in astrophysical experiments. A variety of spectra and 
Susy breaking scenarios have been investigated exper- 
imentally in detail [4]; here just two models, mSUGRA 
and GMSB scenarios, are presented. 



2 t detection & measurement principles 

A typical ILC detector [5] is displayed in Fig. 1. The 
main characteristics, relevant to the present study, are: 
a TPC with excellent tracking and dE/dx resolution 
to identify slow, heavy particles by ionisation; a highly 
segmented hadronic calorimeter (HCAL) with energy 
resolutions 5Eh/E = 0.5/y/E/GcV for hadrons and 
SE em /E = 0.2/y / iJ/GcV for electrons/photons; an 
instrumented iron yoke to allow for muon detection 
and coarse calorimetric measurements of hadrons. The 
amount of material available to absorb a heavy f in the 
HCAL or yoke corresponds to an acceptance for scaled 
momenta of p/m = (3j < 0.4 — 0.5. 
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Fig. 1. Quadrant of a typical ILC detector [5], length units 
in mm; amount of material indicated by R [gem -2 ] 
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The stau detection and measurement principle con- 
sists of several steps: identify a f and determine its 
mass from kinematics; follow the track until it is trap- 
ped inside the detector; observe the stopping point un- 
til a decay f — > tG is triggered by a large energy re- 
lease uncorrelated to beam collisions; record the decay 
time to determine the f lifetime; finally, measure the 
r recoil energy to get the gravitino mass 



The ILC provides a very favourable environment. 
The centre of mass energy can be adjusted to optimise 
the number of observable staus. The e + e _ beams col- 
lide in bunch trains of 1 ms duration repeated every 
200 ms; thus, the detector is most of the time inac- 
tive and in principle ideally suited to measure long- 
lived particles. However, it is envisaged to operate the 
HCAL in a pulsed mode, switching on only during col- 
lisions. Clearly this concept has to be revised to reach 
a reasonable duty cycle. 



3 Experimental analyses - case studies 

The analysis is based on a complete event simulation 
including QED radiation, beamstrahlung and detector 
resolutions. The experimental signature is very clean 
and distinct from Standard Model background which 
can be efficiently rejected. There are no missing par- 
ticles (except v's from decays), the observed particle 
momenta are balanced, | J2iPi\ — 0, but their moduli 
don't sum up to the cms energy J^. pi < y/s. These 
features allow the sparticle masses and decay chains 
to be reconstructed from the event kinematics. Each 
Susy event contains two f's, easily identified by ioni- 
sation in the TPC (dE/dx ~ /3 -2 ), and their passage 
through the detector can be accurately followed. The 
location of stopping f's can be determineded within a 
volume of a few cm 3 . 

The production of low momentum f's with a suit- 
able /?7 factor to be trapped in the detector proceeds 
either directly or via cascade decays from light sleptons 
or neutralinos. All these processes — fifi, e^e^, fj-Rp-R 
and XiXi — rise only slowly above kinematic threshold 
with cross sections a oc /3 3 , thus providing relatively 
low rates. More efficient, if kinematically accessible, 
is associated selectron production e + e~ — ♦ e^e^, in- 
creasing as a (x j3 near threshold. The event signatures 
are multi-lepton topologies: 2fi from pair production, 
2fi2r from ncutralino production and 2fi2r2£ from 
selectron and smuon production. 

3.1 mSUGRA scenario GDM e 

In supergravity mediated symmetry breaking (SUGRA) 
the gravitino mass m 3 / 2 is a free parameter of the same 
order as the other sparticle masses. In minimal versions 
with the f NLSP the common scalar mass mo has to 



Table 1. Sparticle masses and decay modes of mSUGRA 
scenario GDM e accessible at y^s = 500 GeV 
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be small and much lower than the common gaugino 
mass M 1/2 . In the mSUGRA scenario GDM e [6] a 
tighter definition is used with unified scalar and grav- 
itino masses m = m 3 / 2 = 20 GeV, M 1 / 2 = 440 GeV, 
Aq = 25 GeV, tan/3 = 15 and sign/i = +. The corre- 
sponding sparticle spectrum is compiled in Table 1. 

The experimental assumptions for the case study 
are the canonical ILC energy ^/s = 500 GeV and an 
integrated luminosity C — 100 fb~ (less than 1 year of 
data taking). The inclusive f production cross section 
is a(fifiA) = 300 fb. 

The prolific stau production rate is characterised by 
the scaled momentum distribution p/m = shown 
in Fig. 2 a for the various reactions. The majority of 
particles come from diagonal slepton and ncutralino 
pairs and leave the detector (peak around /?7 ~ 1). 
One observes, however, a second peak at low f3^ < 0.5 
from cascade decays of SrCl production, which will be 
stopped in the detector. The number of f's trapped 
are A* cal = 4100 and 7Vf okc = 1850 in the hadron 
calorimeter and yoke, respectively. The choice of en- 
ergy turns out to be optimal. Similar rates can be ob- 
tained at 380 GeV, just above the slepton/neutralino 
thresholds. Selecting fifi pairs to absorb f's is much 
less efficient, e.g. A* cal = 1600 at y/s = 340 GeV for 
the same integrated luminosity. 

The stau mass measurement is based on the kine- 
matics of e + e _ — * fifi, see magenta curve in Fig. 2 a, 
to be identified as a pair of collinear, non-interacting 
particles with momenta pf < yfs/2 = E?. A determi- 
nation of the mean momentum (pf) = 192. 4± 0.2 GeV 
leads to a precise f mass of 

mf = 157.6 ± 0.2 GeV . 

Alternatively one may use the much larger sam- 
ple of all identifed f's, also those leaving the detec- 
tor, and perform a the time-of-flight measurement us- 
ing the calorimeter. The time resolution is St = 1 ns, 
the track lengths arc ~ 2 — 4m. The reconstructed 
mass distribution ttt-toF = 1 P 2 — l)p 2 , displayed 
in Fig. 2 b, provides an accuracy Stutof = 0.15 GeV, 
similar to that of the momentum measurement. 

The stau lifetime measurement is based on the de- 
cays of f's which have been stopped in the detector. 
Requiring an isolated energetic cluster or muon above 
a certain threshold originating somewhere inside the 
sensitive fiducial volume of the calorimeter or yoke, re- 
sults in the decay time distribution shown in Fig. 2 c. 
A fit to the spectrum gives a f lifetime of 

r = (2.6 ±0.05) ■ 10 6 s , 
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mSUGRA scenario GDM e 
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Fig. 2. GDM e scenario, assuming C = 100fb _1 at yfs — 500 GeV: (a) f production spectra of scaled momentum 
p/m = /37 with contributions from various processes; (b) f mass mxoF spectrum; (c) f lifetime distribution; (d) r jet 
energy spectrum of the decay f\ — > tG, compared with simulations of m,Q = 20 GeV, 10 GeV and 30 GeV 



corresponding to roughly one month. For the actual 
conditions the lifetime scales as t ~ (also valid 
for very low masses) and rises significantly faster for 
ttIq > 0.25 m?, reaching 10 8 s for a 75 GeV gravitino. 

Note: The relative precision on the f lifetime does 
not depend on the gravitino mass, should it be much 
lighter as for larger mass splittings or in gauge me- 
diated supersymmetry models. Technically, there may 
be a limitation to measure lifetimes below 1 ms which 
corresponds to the sensitive time of beam collisions. 

A direct gravitino mass measurement can be per- 
formed by exploiting the r recoil of the decay f — > 
tG, see (2). The upper endpoints of the energy spec- 
tra which coincide with the primary r energy E T = 
77.5 GeV, are directly related to the masses involved. 
The leptonic 3-body decays r — > Ivv are not very use- 
ful due to the soft spectrum peaking at low values. Well 
defined upper edges arc provided by the hadronic de- 
cays to heavier final states r — > pv and r — > ttttttv. The 
energy distribution of both decay modes, defined as 'r 
jets', is shown in Fig. 2 d. In order to illustrate the 
sensitivity to the gravitino mass, simulations assum- 
ing the nominal value of = 20 GeV and shifted by 
±10 GeV are shown as well. A fit to the r jet energy 
spectrum, yields a gravitino mass 

m d = 20 ± 4 GeV . 

Combining all results one can test the gravitational 
coupling of the stau to the gravitino and access the 
Planck scale, respectively Newton's constant. Inserting 
the expected values and accuracies on m?, t and m,Q 
in (1) one finds for the supergravity Planck scale 

M P = (2.4 ±0.5) • 10 18 GeV , 

where the error is dominated by the precision on the 
gravitino mass. 

The gravitino mass can be deduced more precisely 
from the f mass and lifetime, if the gravitational cou- 
pling is shown to be responsible for the decay or is 
assumed and the macroscopic value of Mp is taken in 
the decay- width of (1). The resulting gravitino mass 
is mg. = 20 ± 0.2 GeV, where the error is dominated 
by the lifetime measurement. This mass value can be 
used to get access to the supersymmetry breaking scale 



F = V3M P m 3/2 = (8.3 ± 0.1) ■ 10 19 GcV 2 , which is 
an important parameter to unravel the nature of the 
supersymmetry breaking mechanism. 

The expected accuracies on the f and G observ- 
ablcs and derived quantities of the GDM e scenario 
are summarised in Table 2. 

Table 2. Expected accuracies on f and G properties of 
observables and derived quantities; GDM e scenario, based 
on C = 100 fb" 1 at y/s = 500 GeV 

observables derived quantities 
m f 157.6 ±0.2 GeV 

T f (2.6 ± 0.05) 10 6 s rnT 20 ± 0.2 GeV 

m? T 20 ±4.0 GeV M P (2.4 ± 0.5) 10 18 GeV 



It is a unique feature of gravitino LSP scenarios 
that the Planck scale can be directly measured in mi- 
croscopic particle experiments by studying the prop- 
erties of the NLSP and its decay. A further interest- 
ing test to reveal the nature of the gravitino as the 
supersymmetric partner of the graviton would be to 
determine the spin, which is possible by studying cor- 
relations in the radiative decays f — > tjG [7] . 

3.2 GMSB scenario SPS 7 

Gauge mediated symmetry breaking (GMSB) usually 
occurs at rather low scales and a light gravitino is 
naturally the LSP. Typical masses are of order eV to 
keV which may be extended in the GeV range. The 
GMSB reference scenario SPS 7 [8] is described by the 
conventional parameters A — 40 TcV, M m = 80 TeV, 
N m = 3, tan/3 = 15 and sign/i = ±, The sparticles are 
relatively light: m fl = 123.4 GeV, m ijt = 130.9 GeV, 
m; = 262.8 GeV, m^o = 163.7 GeV* The gravitino 
mass is set arbitrarily to nig = 0.1 GeV. 

The SPS 7 model is investigated assuming y/s = 
410 GeV and C = 100 fb~ , with a large inclusive f 
cross section of cr(fifiA) = 420 fb. As seen in the f3j 
distribution of Fig. 3 a, most f's leave the detector. 
There is, however, a large signal at /3-y ~ 0.4 from crcl 
production, contributing to samples of A^ cal = 10000 
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GMSB scenario SPS 7 
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Fig. 3. SPS 7 scenario, assuming C = 100 fb" 1 at yfs = 410 GeV: (a) f production spectra of scaled momentum 
p/m = /?7 with contributions from various processes; (b) f lifetime distribution; (c) r jet energy spectrum of the decay 
fi — > rG compared with simulations of mg = GeV and 10 GeV 



and Nf° c = 4900 trapped f's in the calorimeter and 
yoke, respectively. 

The analysis of fifi pair production yields a mass 
of mfj = 124.3 ± 0.1 GeV. From a fit to the decay time 
distribution, shown in Fig. 3 b, one obtains a lifetime of 
r = 209.3 ± 2.4 s. These values can be used to derive a 
very accurate gravitino mass of ttiq = 100 ± 1 McV as- 
suming a gravitational coupling. To illustrate the sen- 
sitivity to low gravitino masses as expected in many 
GMSB models: a gravitino mass of 0.5 MeV corre- 
sponds a t lifetime of 5 ms, which should be easily 
measurable. 

The t jet recoil energy spectrum is displayed in 
Fig. 3 c. As can be seen from the simulation curves for 
OGeV and 10 GeV gravitinos, the measurement is not 
sensitive to such low masses and can only serve to set 
an upper limit of < 9 GeV (at 95% CL). 

This limitation can be understood from (2). The 
sensitivity to low gravitino masses decreases rapidly, 
e.g. for = 0.1 m? the energy deviates from m^/2 
by one percent, which requires E T to be known with a 
precision well below. Notice that E T is the maximum, 
the upper edge of the observable r jet energy. A direct 
measurement of very large f — G mass splittings be- 
comes experimentally extremely difficult, getting grad- 
ually impossible for ratios m^/mf ^0.1. 

The nature of the LSP remains undetermined with- 
out knowing the gravitino mass. Further information 
can be gained from a study of radiative three-body de- 
cays f — > tjG. The differential decay rates and 7 — t 
correlations for a light spin 3/2 gravitino G to be com- 
pared with a spin 1/2 neutralino x [7] an d a spin 1/2 
axino 5 [9] have been calculated and found to be quite 
different in particular phase space regions of back-to- 
back topologies. Experimentally the analysis of radia- 
tive t decays is quite ambitious: the branching ratios 
are suppressed by two orders of magnitude and single 
7's have to be disentangled from the bulk of photons 
and hadrons in r decays. The excellent performance 
of the proposed high granularity, 'pictorial' calorime- 
ter [5] together with the large data samples expected 
at the ILC should make it possible to discriminate 
between a light gravitino, a neutralino and an axino 
LSP. 



4 Conclusions 

Future ILC experiments have a rich potential to study 
Susy scenarios where the gravitino G is the LSP and 
a charged stau f is the long-lived, mctastablc NLSP. 
Precise determinations of the f mass and lifetime and 
of the G mass appear feasible already with moderate 
integrated luminosity. (More SUSY scenarios can be 
found in [4].) A measurement of the gravitino mass 
from the r recoil spectra of the decay f — > tG gives 
access to the gravitational coupling, i. e. to the Planck 
scale, and thus provides a unique test of supergravity. 
Such observations will put stringent constraints on an 
interpretation of the gravitino as dark matter candi- 
date. 
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